%*******************************************************************************
%*
%*   Copyright ©2009 Alexandru Marin
%*
%*   This file is part of nwave/pMRI C library.
%*
%*   This project started as an MSc degree project during a
%*   stage at the Paris-Est university, under joint supervision of prof.
%*   Bogdan Ionescu (Image Analysis and Processing Laboratory, "Politehica"
%*   University of Bucharest), Caroline Chaux and Jean-Christophe Pesquet
%*   (Laboratoire d'Informatique, Université Paris-Est, Institut Gaspard Monge).
%*
%*   nwave/pMRI C libray is free software: you can redistribute it and/or modify
%*   it under the terms of the GNU General Public License as published by
%*   the Free Software Foundation, either version 2 of the License, or
%*   (at your option) any later version.
%*
%*   nwave/pMRI C libray is distributed in the hope that it will be useful,
%*   but WITHOUT ANY WARRANTY; without even the implied warranty of
%*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%*   GNU General Public License for more details.
%*
%*   You should have received a copy of the GNU General Public License
%*   along with nwave/pMRI C library. If not, see <http://www.gnu.org/licenses/>
%*
%*******************************************************************************



function mri_read_results(filename)

results = mri_read_matrix(filename);

% normalize image
results = 255*abs(results)/max(abs(results(:)));

results_dims = ndims(results);
ncoils = size(results, results_dims);


switch results_dims
    
    case 4 %3d
        
        for i=1:ncoils
           figure;
           image(results(:,:,1,i));
        end
        
    case 5 %3dt
        
    case 3 %2d

        for i=1:ncoils
           figure;
           image(results(:,:,i));
        end

    otherwise
        error('Matrix not recognized');
    
end


end